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ABSTRACT 

The baryonic Tully-Fisher relation is an important observational constraint on cosmological and galactic 
models. However, it is critical to keep in mind that in observations only stars, molecular, and atomic gas are 
counted, while the contribution of the ionized gas is almost universally missed. The ionized gas is, however, 
expected to be present in the gaseous disks of dwarf galaxies simply because they are exposed to the cos- 
mic ionizing background and to the stellar radiation that manages to escape from the central regions of the 
galactic disks into their outer layers. Such an expectation is, indeed, born out both by cosmological numerical 
simulations and by simple analytical models. 

Subject headings: galaxies: kinematics and dynamics - galaxies: dwarf - galaxies: irregular - galaxies: spiral 
-methods: numerical 



1. INTRODUCTION 

The baryonic Tully -Fisher relation (BTFR; iFreemanll 19991: 
iMcGaugh et al.l 120001) extend s the classical Tully-Fisher re- 
lation (ITuUv & Fisher! Il977h by including, in addition to 
stars, all the baryonic content of galaxies. The physical rea- 
son for such an extension is very sound - in galactic halos 
that maintain the universal fraction of baryons, the baryonic 
mass becomes a linear proxy of the total mass. A devia- 
tion from the expected total mass - rotational velocity rela- 
tion may indicate missing baryons (although such deviations 
may be degenerate with the effect of stellar feedback, e.g. 
iDutton & van den Bo sch 2009); hence baryonic Tully-Fisher 
relation becomes a powerful test of galaxy formation models. 

It is not surprising, therefore, that a significant ef- 
fort has been expended over the year s in measuring 
the BTFR in a diverse set of galaxies jVerheiienl 20011: 
Bell & de Jong 2001; Gurovich et al."2004'; McGaugh'2005i 
Pfenniger & Revaz 2005; Begum etal. 2008; Sta rk et al.i 
2009; ' Trachternach et al. 2009; iGurovich etal] 120101: 
Torres-Flores et al.i,201 L). The c urrent state of affairs is well 
svnthesized bv lMcGaughl ( 1201 lb . The remarkable power-law 
behavior of BTFR over 4 decades in mass indicates a 
substantial fraction of "missing" baryons in the lowest mass 
galaxies. The observed "baryonic" components of these 
galaxies are dominated by atomic gas, and the decreasing 
"baryonic" fraction with galaxy mass is often interpreted as a 
substantial mass loss from the dwarf galaxies. 

A limitation of such an argument, however, is that no ob- 
servation actually measures the baryonic contents of galax- 
ies; only stars, molecular, and atomic gas are accounted in 
observational studies, but not the ionized gas (with a handful 
of exceptions). In other words, in dwarf, HI dominated galax- 
ies, the observed BTFR is, in fact, an HI Tully-Fisher relation. 
One needs, therefore, be careful to distinguish the total BTFR 
that accounts for all of the baryons within the virial radius of 
a galaxy, and the observed BTFR that only includes stars and 
molecular and atomic gas, but does not account for the ionized 
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gas. 

In this paper I demonstrate that dwarf galaxies are expected 
to contain warm ionized gas (to be distinguished from hot ha- 
los) - simply because the gaseous disks of galaxies are ex- 
posed to the cosmic background radiation (plus whatever of 
their own ionizing radiation escapes into the outer layers of 
their disks) , which is going to ionize the outer layer of atomic 
disks down the column density A^h 10'^ cm^, comparable to 
the transition column density between Lyman Limit systems 
(which are mostly ionized) and Damped Lyman-a systems 
(which are mostly neutral). 

Hence, the interpretation of the observed BTFR is rather 
non-trivial, nor can it easily be used to deduce the fraction of 
baryons ejected by the feedback from dwarf galaxies. 

2. SIMULATIONS 

The simul ation used in this p a per is similar to the one 
described in iGnedin & Kravtsovl (120101) . Specifical l y, the 
Adaptive Refinement Tree (ART) c ode d Kravtsovl 119991 ; 
iKravtsov et 111120021: iRudd et al.ll2008l) is employed to model 
a 6/z~' Mpc cube centered on a typical 2L, (at z = 0) galaxy. 
The Lagrangian region of a sphere with radius equal to 5/?vir at 
z = (a "region of interest") is sampled in the initial conditions 
with the effective 512^ resolution, while the rest of the simu- 
lation volume is sampled more crudely. This setup results in 
the mass resolution in the region of interest of 1.3 x 10'' 
in dark matter. This region of interest is then allowed to refine 
adaptively as the simulation proceeds in a quasi-Lagrangian 
manner, all the way down to additional 6 levels of refinement 
(total spatial dynamic range of about 30,000), maintaining 
spatial resolution of 260 pc in comoving reference frame. The 
simulation adopts ACDM cosmology similar to the WMAPl 
one (Qm = 0.3, Qb = 0.046, (Ts = 0.9 and h = 0.7). 

The physical processes model ed in the simulation are ex - 
actly the same as described in iGnedin & Kravtsovl (120101) . 
with one exception. In particular, the simulation incorporates 
gas cooling (including cooling on metals, molecular hydro- 
gen, and dust), a phenomenological model for molecular hy- 
drogen formation, full time-dependent and spatially variable 
3D radiative transfer of ionizing and Lyman- Werner band ra- 
diation (both fr om local sources and fro m the incident cosmic 
background of iHaardt & Madauj ( 120011) ) using the Optically 
Thin Variable Eddington Tensor (OTVET) approximation of 
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iGnedin & Abell (l2001h. and star fo rmation in the molecular 
gas usmg iKrumholz & TanI (120071) recipe . The only differ- 
ence from the iGnedin & Kravtsovl ( 120101) simulation is that 
the supernova feedback is disabled. 

Existing models of supernova feedback in cosmological 
simulations (and other types of stellar feedback, except the 
feedback from stellar ionizing radiation, that is explicitly in- 
cluded) are still insufficiently developed to be fully realistic. 
Therefore, I deliberately disable that form of feedback to limit 
numerical treatment to only those physical processes that can 
be modeled with a reasonable degree of realism. 

It is difficult to envision that the supernova feedback would 
increase the baryonic mass of galaxies, hence baryonic frac- 
tions in model galaxies should be considered as upper limits 
on the baryonic fractions in a real universe; any reasonable 
stellar feedback models should reduce these fractions even 
further. 

In order to verify that the results presented below are ro- 
bust and numerically converged, I re-run the last 150 Myr of 
the time evolution in the simulation with an additional level 
of refinement (twice higher spatial resolution) and more ag- 
gressive refinement criteria (to enlarge the regions refined to 
the highest allowed level). This convergence test has demon- 
strated that all the results presented below are highly robust. 

3. TOTAL AND OBSERVED BARYONIC TULLY-FISHER RELATIONS 

In order to compare the simulation to the observations, the 
total baryonic mass (all the way to the virial radius) and the 
observed "baryonic" mass (i.e. the sum of the stellar mass, 
molecular gas mass, and atomic gas mass) are measured for 
all model galaxies from the region of interest. To measure 
the rotational velocity, the mass density profiles for various 
mass components of model galaxies are constructed, and the 
circular velocity Vf are measured at the radius containing 90% 
of all HI. This procedure is similar to the one adopted by 
iMcGaughl (1201 Ih . I have verified that, if the radius including 
either 80% or 95% of all atomic gas is used instead, the results 
would change by no more than a symbol size in Figure [T] 

Since the region of interest is centered on a single large 
galaxy, only a handful of other noticeable isolated galaxies 
end up in that region by z = 0. To somewhat enlarge the sam- 
ple, I show in Figure[T]both z = and z = 0.5 simulation snap- 
shots, which are apparently fully consistent with each other 

Figure [T] shows the two BTFR (the total baryonic one and 
the observed BTFR that only includes stellar mass, molecular 
gas mass, and atomic gas mass but not the ionized gas mass). 
The s olid line tracks the best-fit observed BTFR (iMcGaughl 
1201 lb relation in the form 
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The dotted line gives a relation between the product of the uni- 
versal baryon fraction funi = fls/^M and M200 (mass within 
the overdensity 200 with respect to the critical density) vs 
Vmax for dark matter halos f rom Bolshoi and Multi-Dark sim- 
ulations (iPrada et al.ll201 ll) . This relation is slightly steeper 
than /uni^aoo oc V^^^ '^^^ 1° 'he halo concentration depen- 
dence on halo mass. The simulation lies to the right of that 
line in Figure [T] because Vmax from pure N-body simulations 
is not a sufficiently accurate proxy for the disk circular veloc- 
ity - halo contraction due to baryonic condensation results in 

y/«(i.2-i.3)ymax. 

It is also possible to estimate the HI masses of dark matter 
halos using abundance matching - i.e. identifying dark mat- 




V, (km/s) 

Fig. I. — The baryonic Tully-Fisher relation for model galaxies at z = 
(black) and z = 0.5 (gray); I show two redshifts since the z = sample is very 
small by construction. Each model galaxy is represented by a vertical line 
with two symbols: open square shows the total baryonic mass of the galaxy; 
filled square shows the observed "baryonic" mass, i.e. the sum of stellar mass, 
inolecular gas inass, and atoinic gas mass. The difference between the two 
squares is the contribution of the ioniz ed gas. The solid line is the best-fit 
observed BTFR froin [McGaughl 1201 IH ; the dotted line is /uniM200 vs V^ax 
for dark matter halos in the ACDM cosmological model, while the dashed 
line is Mm vs Vmax obtained by abundance matching of the theoretical halo 
mass function and the observed HI mass function (see text for details). 



ter halos from a theoretically computed halo mass function 
with HI galaxies fr om the actual observed HI mass function 
dZwaan et al.ll2005h that have the same spatial number den- 
sity. The details on how thi s matching is performed are de- 
scribed in lMarm et al] ( 120101) or Trujillo-Gomez et al. (201 1); 
HI mass vs Vmax that results from such abundance matching 
is shown in Fig. [T]with the dashed line. Not surprisingly, it 
matches the best-fit observed BTFR at low masses very well. 

The observed BTFR in the simulation (solid squares) is not 
far from the best-fit line to the observational data, but is above 
it by about a factor of 2 for smallest galaxies. It is important 
to remember that the simulation should only be treated as an 
upper limit to the baryonic content of galaxies, because the 
stellar feedback is deliberately disabled in the simulation (ex- 
cept the feedback of ionizing radiation, which is included in 
the simulation by virtue of following the full spatially-variable 
and time-dependent 3D radiative transfer). The fact that the 
simulation points are not falling on the observed relation for 
Vf < 200 km/s simply illustrates the well-known fact that the 
stellar feedback processes are important for determining the 
baryonic contents of sub-L* galaxies. 

In order to illustrate the crucial difference between the to- 
tal baryonic fraction and the observed (stars -1- molecular -1- 
atomic gas) baryonic fraction in model galaxies, I show the 
two fractions in Figure |2l As one can see, the total baryonic 
fraction in model galaxies remains close to universal all the 
way to Vf ^ 40 km/s, and at lower Vf the effect of photo- 
evaporation of gas from dark matter halos due to heating by 
the cosmic ionizing background radiation becomes important. 
That effect in thi s simulat i on is approximately described by a 
fitting formula of lGnedinI ( l2000h . 
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with parameter values a = I and Mg = 7 x lO^MfT). Th e lat- 
ter value is the same as found bv lOkamoto et al.l ( 120081 con- 
verted from the virial overdensity of 97 with respect to the 
critical used there to the overdensity of 200 used here, for 
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Fig. 2. — The baryonic fraction (in units of the universal fraction = 
^b/^m) for model galaxies at j = (black) and z = 0.5 (gray). Symbols are 
the same as in Figure [T] The dotted line approximately traces the effect of 
photo-evaporation of gas from halos due to cosmic background (Equation [T) 
while the short-dashed line shows a simple model (3) for the neutral content 
of galactic disks in dwaif galaxies. The long-dashed line is the fraction of 
ionized gas within the HI radius (Equation|4). 

the halo concentration of 14)o This function (plotted against 
Vmax) is shown in Figurep as a dotted line. In agreement with 
iMcGaugh & Wolll (120101) . this line alone is unable to explain 
the baryonic contents of dwarf galaxies. 

It is also possible to construct a simple analytical model 
that approximately reproduc es these results f or the observed 
BTFR, following the ideas of lMo et al] (119981) . Let's consider 
an exponential gaseous disk with the surface density run 

such that the disk scale length is a given fraction /i of the halo 
virial radius R2Q0, 

Rd = fiRioo, 

and the disk mass is a fraction ly of the total halo mass of M200, 
so that 

(expressions for /i and ly are given in lMo et aI1(ll998h as func- 
tions of other parameters). The disk is exposed to the cos- 
mic ionizing background, which ionizes hydrogen (from both 
sides) in a slab of gas below the surface density Ehii- Then 
the mass of the neutral gas in the disk is 

/•^HII 

MHi = 27r/ (j:(R)-j:Hu)RdR, (2) 
Jo 

where Run is the edge of the HI disk (I](7?hii) = Shu)- The 
minus Shu term under the integral in Equation (|2|i appears 
because both sides of the disk are exposed to the ionizing ra- 
diation, and hence each side has an ionized layer of thickness 
Ehii/2. Integral |2] can be easily taken, so that the fraction of 
the disk that is neutral is 

VM200 So y Shii 2 V Shu/ J 

* I also notice that several recent studies of the photoioniz ation effect 
that did not treat the radiative transfer of ionizing radiation IHoeft et"al] 
12003 : ICTrnneTan iTOO?; Okamoto et al. 2008) found that the best-fit value 
of a = 2. However, in the simulations with the full spatially-variable and 
time-dependent 3D radiative transfer we always find that the best-fit value of 
o is 1 , even in simulations run with widely different codes, like the moving 
mesh code of [Gnedin ( 1 9951) and the currently used ART code. 
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Fig. 3. — The circular and rotational velocity profiles (top) and surface 
density profiles for the total, ionized, and neutral gas (bottom) for two model 
galaxies at z = 0, one with Vf = 40 km/s (black lines) and another one with 
Vf = 23 km/s (gray lines; the lowest of two black symbols at this velocity in 
Figures[T]and|2). In the circular velocity plot, the solid segment of the circular 
velocity shows its extent over the H I disk, whereas the dotted line shows the 
extension of the line into the region where there is no atomic gas. 

This dependence (for the combination ly/fJ? = 8, which is 
somewhat on a lower side but still plausible, and Shu = 
0.4Mo/pc2 - the reason for the latter choice will become 
clear belowQ), multiplied by Equation ([T]) and plotted as a 
function of Vmax, is shown in Figure |2] with the dashed line. 
It does not fit the simulation results perfectly, but roughly 
traces the lower envelope of filled symbols; this simplistic 
model does not account for stellar contribution to the observed 
BTFR, and is therefore can only serve as an approximation to 
the lower limit of the observed baryonic fraction of model 
galaxies. 

From an observational perspective, it may be useful to know 
what fraction of the ionized gas is actually inside the HI disk. 
In the simple model this quantity is easy to compute, 

Mmi(R<Rmi) = 2T: ^miRdR, 
Jo 

or 

. , MmiiR < Rmi) IShh/, V 

fmi(R < Rmi) = = In-^ • (4) 

iyM2Qo 2 So V Shu/ 

The fraction of ionized gas outside the HI disk is then simply 
one minus the sum of Equations (O and Equation (|4]i. The 
long-dashed line in Figure|2]shows Equation (HJi as a function 
of Vmax- The ionized gas starts dominating the disk mass only 
for Vmax < 30 km/s, but makes a > 20% contribution all the 
way to Vmax ~ 100 km/s. 

Finally, one may ask where the ionized gas, apparent in Fig- 
ures [1] and |2] is actually located in model galaxies? To an- 
swer that question, I show in Figure [3] surface density profiles 
for the total, ionized, and neutral (atomic and molecular, al- 
though the molecular fraction in the shown galaxies is small) 
gas in two galaxies, in which the ionized gas contribution is 
large. In both galaxies neutral gas forms an approximately ex- 
ponential disk, while the surface density of the ionized com- 

^ This is also the column density of a typical sub-DLA absorption system. 
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ponent remains approximately constant well outside the scale 
length of the neutral disk. These trends are consistent with 
the simple analytical model presented above. The choice of 
Shii = O.4M0 / pc~ is the actual value of the fixed surface den- 
sity of the ionized gas found in the simulation. 

The top panel of Fig.[3]also shows the circular velocity pro- 
files for the two galaxies and the rotational velocity of H I gas. 
The difference between the two is due to non-circular mo- 
tions in the gas; t his difference is consistent with the prior 
theor etical models ( Valenzuela et al.l l2007l : lDalcanton & Stilpl 
120 lot FTruji llo-Gomez et a l. 2()l_y and observational mea- 
surements of galactic rot ation curves (c.f. iTrachternach et alJ 
IM)8tlAdamsetaDl20Tlh . 

4. CONCLUSIONS 

The baryonic Tully-Fisher relation is an important observa- 
tional constraint on cosmological and galactic models. How- 
ever, it is critical to keep in mind that in observations only 
stars, molecular, and atomic gas are counted, while the contri- 
bution of the ionized gas is almost universally missed. Hence 
the observed BTFR does not count all baryons. Comparison 
of such observations to theoretical predictions is, therefore, 
highly non-trivial, and requires a proper modeling of radia- 
tive transfer of ionizing radiation, at least in an approximate 
form. 

In this paper I present an example of such modeling in the 
form of a cosmological numerical simulation with radiative 
transfer The simulation is not fully realistic, since it does not 
include a treatment of supernova feedback, and thus overesti- 
mates the baryonic fraction of model galaxies. Even with this 
incomplete treatment, low mass model galaxies contain large, 
and for Vj < 50km/s dominant, contribution of ionized gas. 

The ionized gas is present in the gaseous disks of model 
galaxies in two regions: the outer parts of the disks are ion- 



ized since their surface densities (and, hence, column den- 
sities) are too low. Even more importantly, the outer layers 
of the inner disks are also ionized, because they are exposed 
to the cosmic ionizing background and to the stellar radiation 
that manages to escape from the central regions of the galactic 
disks into their outer layers. These layers are direct analogs 
of the Reynolds ' layer observed in the Milky Way galaxy 
(lRevnoldslll993l) . but their relative contribution to the total 
mass budget becomes progressively larger as galactic disks 
become less massive, less dense, and allow ionizing radiation 
to reach deeper. 

The existence of this ionized gas is not just a theoretical 
conjecture - it is unavoidable from purely physical grounds, 
since the cosmic ionizing background must ionize the outer 
layers of galactic HI disks (on both sides) down to column 
densities of Lyman Limit systems, A^hi ^ 10'"^ cm"- (in my 
simulation it is Shii = O.2M0/pc^ or, equivalently, HI = 
2.5 X 10''*cm-2). 
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